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Abstract 

We give a dynamical justification of a differential cross section for the Schrodinger 
equation in the context of long time transition to stationary regime. Our approach is 
based on spherical incident waves, produced by a harmonic source, and uniform long-range 
asymptotics for the corresponding spherical limiting amplitudes. The main result is the 
convergence of the spherical limiting amplitudes to a limit as the source is moving to infinity. 
The main technical ingredients are the Agmon-Jensen-Kato analytical theory of the Green 
function and the Ikebe uniqueness theorem for the Lippmann-Schwinger equation. 
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1 Introduction 

The differential cross section is the main observable in quantum scattering experiments. The 
notion was introduced first to describe the Rayleigh scattering of sunlight and the Rutherford 
scattering of alpha particles as the quotient 

°(d)=fa c (d)/\r\ • (i.i) 

Here j m is an incident stationary flux, and jl c (9) is the angular density of the scattered stationary 
flux j sc (x) in the direction 9 G M 3 , \9\ = 1: 

3l c (9)= km R 2 f c (R9)-9 , (1.2) 

R— >oo 

where j sc (x) denotes the scattered stationary flux (see Fig. [T]). 









• 







Figure 1: Incident flux and scatterer. 

In both scattering processes, the differential cross section is well-established in the framework 
of the corresponding dynamical equations: the Maxwell equations in the case of the Rayleigh 
scattering and the Newton equations in the case of the Rutherford scattering. On the other 
hand, a satisfactory justification of the quantum scattering cross section is missing. 

Our goal is to complete the theory of the quantum differential cross section in the framework 
of the Schrodinger equation, 

iip(x, t) = Hif>(x, t) := ~Ai/;(x, t) + V(x)ifi(x, t) , x G M 3 . (1.3) 

The corresponding charge and flux densities are defined as 

p{x,t) = \^j{x,t)\ 2 , j(x,t) = Jm[i/>{x,t)Vif){x,t)] . (1.4) 
The densities satisfy the charge continuity equation (see [23J) 

p(x,t) +div j(x,t) = 0, (i,i)6l 4 . (1.5) 
We justify the formula for the differential cross section 

a{k,9) = 16ir A \T{\k\9,k)\ 2 , (1.6) 



which is universally recognized in physical and mathematical literature: see, for example, 
ESI ED I2H EZ] • Here k G M 3 is the "wave vector" of the incident plane wave e tkx and the T -matrix 
is given by 

T(k', k) := j^(T(E k + i0)e ikx , e lk ' x ) , k', k G M 3 , (1.7) 
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which is the integral kernel of the operator T(E k + iO) := V — VR(E k + iO)V in the Fourier 
transform. It is well known that the T-matrix allows us to express the integral kernel S(k', k) 
of the scattering operator S in the Fourier transform, see [3l [191 121] : 



S(k', k) = 5(k' - k) - m8{E k , - E k )T(k', k), k', k e 



'U 



The problem of dynamical justification of the differential cross section is suggested by the 
discussion in [2TJ pp. 355-357]. Namely, the commonly used "naive scattering theory" consists 
of the following statements 



I. The incident wave is identified with the plane wave 



ip' m (x,t) = e i{kx - Ekt) 



E k = -k 



k e 



;i.9) 



The wave propagates in the direction of the wave vector k and is a solution to the free Schrodinger 
equation (jl.3p with V(x) = 0. 
















»■ 

\j/ in 









Figure 2: Incident plane wave and outgoing spherical wave. 
II. The corresponding solution "intuitively" is identified as the long time asymptotics 

ip(x,t) ~ A(x)e~ iEkt , t^oo, (1.10) 
where the amplitude A(x) is expressed by 



A(x) = e lkx - R(E k + iO)[V(x)e 



\kx~\ 



[l.ll) 



This amplitude admits the following "spherical" long range asymptotics [H p. 279, formula 
(3.58)] 

p *|fc|-|x| 

A(x) ^ e lkx + a{k,e)—— , |x|^oo, 9 := x/\x\ ; (1.12) 



x 



see Fig. [21 
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III. By (H3|), asymptotics f lTT2|) give 

j in = fc, jr(0) = |«(M)| a |*l , (1.13) 

and hence the differential cross section reads 

a(k,9) = \a{k,6)\ 2 . (1.14) 
It is well-known that a(k,9) is proportional to the T-matrix [211 formula (97a)], 

a(k,9) = -4n 2 T{\k\e,k) . (1.15) 

Hence, ([LTD reads as fOD . 

A heuristic derivation of relations (11.101) . (11.111) can be found in [2"Tl pp. 355-357]. However, 
a mathematically consistent treatment of the relations ( in a time dependent picture was not 
proved until now. Moreover, relation (11.61) is considered as a definition of the differential cross 
section: see [EH formulas (1.2), (A.1.6)], [2H formula (96)]), and [281 Definition 7.9 on p. 254]. 

The main problem is mathematical justification of (11 .101) and (II. lip is related to the lack of 
a mathematical model for incident wave tp m (x,t) which would provide convergence (ll.lOp to a 
stationary regime, and at the same time satisfy the "adiabatic condition" 

ip' m {x,t) ->• 0, t-> -oo, xeM 3 , (1.16) 

which corresponds to the spirit of scattering theory. The plane incident wave ( 11. 9ft in the "naive 
scattering theory" does not satisfy ( 11.161) since the wave occupies the entire space. The plane 
wave is a solution to the free Schrodinger equation 

iip(x,t) = ~Ai[)(x,t) , xgK 3 . (1.17) 

Adiabatic condition (I1.16P in acoustic scattering is provided by the "semi-infinite" incident plane 

wave 

</> in (x,t) = Q(\k\t - kx)e i{kx ~ m 

for t < 0, where is the Heaviside function. The incident wave is a solution to the acoustic 
equation 

rb(x,t) = Aip(x,t) , xER 3 \B, (1.18) 

for t < —D if the scatterer is located in the region \x\ < D. The similar incident plane wave can 
be constructed for the Maxwell equations, which makes apparent the meaning of the differential 
cross section in the Rayleigh scattering. 

On the other hand, a similar semi-infinite incident plane wave does not exist in the case of 
the Schrodinger equation. Indeed, we could fix R ^> \k\D and take the semi- infinite plane wave 

vh in (x) = G(-R - kx)e ikx 

as the initial condition at t — 0. However, the corresponding solution does not satisfy the 
adiabatic condition for t — > — oo. The problem is of the great importance also in the context of 
quantum field theory, where the incident and outgoing plane waves play the fundamental role 
[El |22l [231127]. 
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In the traditional approach, the incident wave is a specific initial field which is a solution to 
the corresponding free wave equation in the entire space. On the other hand, in practice, the 
incident wave is a beam of particles or light produced by a macroscopic source and satisfies the 
free wave equation only outside the source. One could expect that, for a large time, the incident 
wave near the scatterer will be an asymptotically free plane wave if the source is "monochro- 
matic" and distant from the scatterer. This model corresponds obviously to spherical incident 
waves, which are standard devices in optical and acoustic scattering [3]. We justify relations 
f ri~T0|) . dUIIJ, and jnSP, dHH in the following steps: 

A. First, we prove a limiting amplitude principle for the harmonic source, i.e., long time con- 
vergence to a stationary harmonic regime with a "spherical limit amplitude", which does not 
depend on initial state. 

B. Second, we prove the convergence of the spherical limit amplitudes to a plane limit amplitude 
when the source goes to infinity: D — > oo. 

Convergences A and B imply that relations ( ll.lOp . ( II. lip , and HI . 13 j) . dl . 14[) hold true in this 
double limit: first t — > oo and then D — > oo. 

Let us comment on our techniques. We use Agmon-Jensen-Kato theory of the resolvent (or 
the Green function) [JJ HH] to prove the limiting amplitude principle, uniform bounds and long 
range asymptotics for the spherical limiting amplitudes. The proof of convergence B relies on 
an application of the Ikebe Uniqueness Theorem for the Lippmann-Schwinger equation. We 
refine the known estimates [31 Lemma 3.2 in Chapter 4] for the derivatives of the remainder in 
asymptotics HI . 1 2j) . 

Note that the flux in our approach is given by the second formula (jl.4j) . It describes the wave 
flux in the one-particle Schrodinger equation ( II. 3p . Respectively, a many-particle interpretation 
of cross section (11. ip is not straightforward. An alternative approach is based on a stationary 
random process for an infinite number of particles modeled as normalized wave packets [8] ITO"! 126] . 
It is instructive to note that both approaches give the same expression for the differential cross 
section. One could expect that a unified approach, combining both a stationary nonrandom flux 
with a many-particle interpretation, should be developed for the second quantized models [5]. 

Let us comment on the known arguments for formula ( II .6p . In the physical literature the 
formula is usually justified by squaring the Dirac delta function and the subsequent physical 
interpretation of arising divergence; see, e.g., [HI [22j [231 12Z] • Another physical approach [23] is 
based on random incident wave packets ip m (x, 0) which are asymptotically proportional to the 
plane waves e lkx : 

\4> in (k', 0)| 2 -» 8{k' - k) . (1.19) 

The known mathematical justifications rely on Dollard's fundamental result [7] on scattering 
into cones. For example, Dollard's result is used in [26] for a clarifying treatment of formula 
( II. 6p . Namely, the normalized angular distribution of a finite charge, scattered for infinite time, 
converges to the normalized function (II. 6p in limit (I1.19p . 

There are also mathematical justifications of the approach [24J. Namely, Dollard's result 
was refined in [61 [15] and [21 Section 3-3] , where the flux across the surface theorem is proved. 
This result was later developed in j8j [9J [Til 125] and applied for a justification of formula (I1.6P in 
context of the Bohmian particle mechanics and incident stationary random processes constructed 
of normalized wave packets (I1.6P in limit ( I1.19p . The survey can be found in [TO] . 
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Let us stress that previous results do not concern the long time transition to a stationary- 
regime. Moreover, we do not exclude the discrete spectrum of the Schrodinger operator H in 
contrast to |8j. 

Our paper is organized as follows. In Section [2] we state the main results. In Section [3] we 
establish the limiting amplitude principle. In Section H] we obtain long range asymptotics and 
uniform bounds for the limit spherical amplitudes. In the next two sections we prove convergence 
B and the corresponding convergence for the flux. Finally, in Section [S] we check formulas (11 . 15)) 
and ffl~T3"j) . which justifies ( Q3|) and fTTB]) . 

2 Main results 

We consider the Schrodinger equation with the harmonic source: 

iip(x, t) = Hip(x, t) - p q (x) e~ iEkt , t > 



x EM. . (2.1) 



Here H = — |A + V(x), = k 2 /2 for fc £ I 3 \ 0, and p q (x) := — q) is the form factor of 

the source. The factor \q\ is introduced for a suitable normalization. 

Let us define weighted Agmon-Sobolev spaces T-L s a = Ti^iM. 3 ) with s, o 6 R. 

Definition 2.1. T-L s a = ?^(1R 3 ) denotes the Hilbert space of tempered distributions ip(x) with the 
finite norm 

U :=||(V) S ^)|U <oo. (2.2) 



We will assume the following conditions. 
HO. The initial state is a function from the space %° ao with some cr > 5/2. 
HI. For some e' > 

sup(a;) 4+£ '|p(a;)| < oo . (2.3) 

2 SIR 3 

H2. The following Wiener condition holds: 

p(\k\d) := J e iW0x p(x)dx ^ , E M 3 , \9\ = 1. (2.4) 

H3. The potential V(x) is a real C 2 function satisfying the condition 

sup(x) 5+£ \d a V(x)\ < oo , M<2, (2.5) 



with some e > 0. 

Finally, we introduce our key spectral assumption. Denote 

M a := {i/j e C\ ■ tfj + Ro(0)Vijj = 0} , 

where -Ro(^) : = (~ |A — u)^ 1 is the free resolvent. The space M. a = M. does not depend on 
o G (1/2, (5 + e)/2) by the arguments before Lemma 3.1 in [16J. 
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H4. We assume: 

Spectral Condition : M = . (2.6) 

The condition holds generically, see the discussion before Lemma 3.1 in 



Let us outline our plan. 

I. First we will prove the limiting amplitude principle: 

N 

^(x, t) ~ <p q (x, t) = B q {x)e- iEkt + C^xy-^ 1 , t^oo, (2.7) 

i 

where ipi{x) denote the eigenf unctions of H corresponding to the eigenvalues ui < 0. The 
asymptotics hold in 'H 2 _ r7 with a > 5/2, and the limit amplitude B q (x) is given by 

B q (x) = R(E k + tO) Pq , (2.8) 

where R(u) := (H — uj)^ 1 is the resolvent of the Schrodinger operator H. 

The coefficients C\ depend on the intitial state ip(x, 0). On the other hand, it is crucially 
important that the coefficients C\ do not depend on q, while the eigenfunctions ipi{x) rapidly 
decay at infinity by Agmon Theorem [U Theorem 3.3]. Hence, the sum over the discrete spectrum 
on the RHS of (12. 7p does not contribute to the scattering cross section. This is why we will omit 
the sum in all calculations below for the simplicity of notation (e.g., as in the case of empty 
discrete spectrum). 

II. Second, we denote B D (x) := B qD (x), where qo = —nD with n = k/\k\ and D > 0, and 
establish the following "spherical version" of long range asymptotics (j!.12p : 

Itol j\k\(\a>- qD H gD \) + aD ( kj9 f 



B D (x) ~ b D (n) 

\x — qu\ ' ' ' \x\ 

as \x — qo\ — > oo , \x\ — > oo , (2.9) 
where 9 := x/\x\ and &d(ti) := 6(n)e l ' fc ' D with b(n) ^ 0; see Fig. [31 

III. Further, we prove the long range convergence of the spherical limit amplitudes, which is our 
central result: for k ^ 

A D {x) := B D {x)/b D {n) ->• A(x) , D -> oo , (2.10) 

where A(x) is expressed by (11. lip . 

IV. At last, (12.101) implies the convergence of the corresponding limit solutions ip£>{%,t) '■— 

(p D (x,t)/b D (n) ->• A(x)e~ lEkt , D oo, (i,i)eK 4 , (2.11) 

and of the corresponding flux (11.41) : 



j D {x) = lm[(p D (x,t)V<p D (x,t)] — >• j^x) = \b(n)\ 2 lm[A{x)VA(x)} (2.12) 

D -> oo 
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Figure 3: Incident and outgoing spherical waves. 



V. Finally, we show that convergence (I2.12p and formula (11. lip justify (ll.ip . ( II. 6p in the limit 
D -)■ oo. 

We establish long range asymptotics (12. 9p using a "spherical version" of the Lippmann- 
Schwinger equation and applying Lemma 3.2 from [3] together with the properties of the resol- 
vent of the Schrodinger equation obtained in [16]. The asymptotics imply convergence (I2.10p 
by a compactness argument and the Ikebe Theorem on the uniqueness of the solution to the 
Lippmann-Schwinger equation [3j [H] which holds under the condition k ^ 0. One of the key 
observations is that the spherical incident wave from (12 .9p becomes the plane incident wave from 
(I1.12p asymptotically when the source goes to infinity: 

, ^ j\k\(\x-q D \-\q D \) _^ e ikx^ D^OO . (2.13) 

F - Qd\ 

Remark 2.2. We suggest that our model of the incident wave suits the physics of quantum 
scattering. Namely, an incident electron beam is produced by a time periodic source like a 
heated cathode in an electron gun. The cathode is described by a density factor, and the source 
is not at infinity though its distance from the scatterer is sufficiently large. 



3 Limiting amplitude principle 



We deduce the limiting amplitude principle (12. 7p from the dispersion decay in weighted energy 
norms established in |16j . 
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Lemma 3.1. Let conditions H0-H4 hold and k G M 3 . Then the limiting amplitude principle 
(12. 7p holds in the norm ofW}_ a with any a > 5/2, where the coefficients Ci do not depend on q 
and the limit amplitude is given by (12. 8p . 

Proof. We should prove that 

N 

il>(x,t) = B q (x)e- iE ^ + J2 C M^~ Mt + r(x,t) , (3.1) 

i 

where 

ll r (")^)lkl -> j t -> oo . 
The solution to the Cauchy problem (12.11) is unique and given by the Duhamel representation 

^(t) = U(t)ip - i [ e~ iEkS U(t - s)p q ds . (3.2) 
Jo 

Here 

N 

U{t)ip Q = Ci^ie' luJlt + r (t) , (3.3) 
i 

where Ci do not depend on q and 

\\r (t)\\ H ^<C(ty^ (3.4) 

by the dispersion decay [H3, (10.9)], which holds in our case due to HO and H3-H4. On the other 
hand, 

I(t) := [ e~ iEkS U{t - s)p q ds = e~ iEkt [ e iEkT U(r) Pq dT , (3.5) 
Jo Jo 

where p q G H®, with some o' > 5/2 by HI. Hence, 

\\U(r) Pq \\ Hla <C(r}-^ (3.6) 

again by the dispersion decay [16, (10.9)]. Therefore, 

I(t) ~ B q (x)e- iEkt , t -»• oo , (3.7) 

where the asymptotics hold in % 2 _ a ■ In particular, contribution of 1(f) into the coefficients C; 
vanishes. Now (13.11) follows from (j3.2p - (13.7p with the limit amplitude 

poo /*oo 

B g (x) = -i e lEkT U(r) Pq dr = -i e i{Ek+i0)T U(T)p q dr . (3.8) 
Jo Jo 

Finally, this integral can be written as (12. 8p . □ 
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4 Spherical waves 

Here we obtain long range asymptotics (|2.9p . Denote R = R(E k + iO) and R = Ro(E k + iO), 

where Ro(uj) = (H — uj)~ 1 is the resolvent of the free Schrodinger operator H = — —A. Let 

us rewrite formula (12.81) for the limit amplitude as the following "spherical version" of the 
Lippmann-Schwinger equation: 

B q (x) = R oPq (x) - R VB q (x) (4.1) 

since R = Rq — RqVR. The free Schrodinger resolvent Rq{u) = (Hq — uj)^ 1 is an integral operator 
with the kernel 

e i>/2u\x—y\ 

R {u,x,y) = — r, weC\[0,oo). 

2tt\x — y\ 

Hence, Rq is the integral operator with the kernel 

e i\k\\x-y\ 

R (E k + iO, x, y) = — , (4.2) 

2tt\x — y\ 

since y/2(E k + iO) = \k\. 

For the first term on the right-hand side of (14.11) asymptotics (I2.9P follows by Lemma 3.2 
from J3j Chapter 4]: 

Lemma 4.1. Under condition HI the asymptotics 

R p t x ) = b( ^ e m\-i\ + K ( x _ q \ xgR 3 , (4.3) 

V \x — q\ J 1 + \x — q\ 

hold, where 

b{6) = V2^p(\k\6) , |0| = 1 , (4.4) 

\K(x - q)\ < C\q\(l + \x- ql)' 1 ^' . (4.5) 
As a corollary, we obtain the bound 

\R p q (x)\ < : | , x e R 3 . (4.6) 

1 + \x — q\ 

For the second term on the right-hand side of (14. ip we need two additional technical lemmas. 
Lemma 4.2. Under conditions HI and H3 the following bound holds: 

sup || V\B g (:r)||£2 < oo for any a < 2.5 + s . (4.7) 

q£R 3 

Proof. The Lippmann-Schwinger equation (14. ip implies 

(1 + VR )VB q = -VR oPq . (4.8) 
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On the other hand, (1 + VRq)- 1 = 1 - VR. Hence, 

VB q = -(1 - VR)VR oPq = -VR oPq + VRVR oPq 

Let us estimate each term on the right-hand side separately, 
i) Condition (I2.5P with a = and bound (14.61) imply 

\VRo Pq (x)\ < — 



(4.9) 



C\q 




(1 + 


x-q\)(l + 


x|) 5 + £ 



x G 



(4.10) 



Therefore, 



\VRo Pg (x)\ < 



C 



x G 



(4.11) 
(4.12) 



(1 + |x|) 4 + £ 

Hence, 

VR p q eCl , a < 2.5 + e . 

Thus the bound (14. 71) holds for the first term on the right-hand side of (14.91) . 

ii) It remains to estimate the last term of (14. 9ft . By (14. 12ft we have RVRop q G ^ 2 CT1 for any 
<j\ > 1/2 since the resolvent R = R(Ek + i0) : C? — > C?L ax is continuous by [161 Theorem 9.2] 
because E k > 0. Therefore, VRVR p q G C\ for a < 4.5 + e by ([23} with a = 0. □ 



Lemma 4.3. Under conditions HI and H3 t/ie following uniform decay holds: 



sup \RqVBJx)\ < c(i + \x\y 



x G 



Proof. We adopt the arguments from the proof of Lemma 3.2 in |3, p. 275]. Namely, 



R VB q (x)\ < C 



\VB q {y)\ 
\x - y\ 



dy 



(4.13) 



C 



\VB q (r,v,0)\dp sin Odd 



r 2 dr 



\J\x\ 2 + r 2 — 2\x\r cos 6- 
by (14. 2p . Applying the Cauchy-Schwarz inequality to the inner integral, we obtain 



\RoVBJx) 



< C 



\VB q (r,^,9)\ 2 d^sm9d9 



d(p sin 6d9 



Ixl 2 + r 2 — 2\x\r cos 8 



r dr 



(4.14) 



C 



\VB q (r,(p,9)\ 2 d(psm6d6 



\x\r 



log 



\x 


_l_ r 


\x 


— r 



r 2 dr. 



(4.15) 
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Applying the same inequality to the last integral, we obtain 

\RoVB q (x)\ 



< C 



x 



(1 + r) 2a \VB q (r, tp, 6)\ 2 d<p sin 6d0r 2 dr 

\t\ 4- r r 2 

log 



JO 



x 


_l_ r 


\x 


— r 



\2u 



<C(a) 
= C(a) 



log 



log 



X 


_l_ r 


\x 


— r 



\2cr 



1 + s |x|sds 



, |l-s| (1 + s|x|) 2ct 
by the uniform bound (14.71) . Splitting the region of integration 

(0,oo) = (0,~)U(~ ~)U(~,oo) 

we estimate the last integral by 

* 1/2 \x\s 2 ds 



Ci{a 
= CAc 



o (1 + sN) 



2<T 



+ (l + |x|) w<7 + 





a; (is 


(1 + 5 





"M/2 r 2j 

2 ' +(l + |a;|) 1 - 2CT + 



\2<T 



31x1/2 



(1 + r) 2CT . 



, (1 + r 
<C 2 (a)(l + |a;|) 1 - 2ff 
if 3 — 2cr < 0. Hence, (I4.13P follows since we can take any a < 2.5 + e by (14. 7p . 

Now we are ready to prove (12. 9j) . 
Proposition 4.4. Asymptotics ( j2. 9f) /ioW under conditions H1-H3. 
Proof. The Lippmann-Schwinger equation (14. ip yields 

VB,(x) = -W2o/> ff (x) - VRoVB q {x) . 
Hence, ([23]) with a = and fICTD . fT4~T3"]) imply 

Therefore, similarly to (14. 3ft . we obtain the asymptotics 

i2o^5,(x) = c ( — )^— r-r + L (x 



X / 1 + X 
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where 

\L q (x)\ ^Ca + M)- 1 - 5 . (4.18) 

Now flHD and (14371) imply 

B,(«) ~lf^ JiL^W +ef f * ^ (4.19) 



x — 5| / \x — q\ \\x\/ \x 

as \x — g| — > oo and |x| — >■ oo. Denote Bd{x) := B qD (x), where 5d = — nZ) with n = k/\k\ and 
L> > 0. Then 

B D (x) ~ 6(n) f , . e^ 11 ^ 1 + dnik, 0)^-^-1 (4.20) 

L|x — 5d| |x| J 

as |x — 5d| — >■ oo and |x| — >■ oo, where 6* := x/\x\ since o(n) ~ p(\k\n) ^ by (14. 4p and the Wiener 
condition H2. This is the unique point in our proofs, where the Wiener condition is needed. 
Finally, (jOOj) can be written as (Q with o D (n) := 6(n)e i l*l £ ' and ar>(Jfe, 0) = #)e" i|fc|D . □ 



Corollary 4.5. Bound (14 ,6p . asymptotics (14.17j) . and formula (14. ip imply 

|B g (g)l < : I ~l~ / . ~i i n , x,gGM 3 . (4.21) 



5 Plane wave limit 

Here we obtain convergence (I2.10p from the uniqueness of solution to the Lippmann-Schwinger 
equation 

A(x) = e ikx - R VA(x) , (5.1) 
which is equivalent to (11.1 ip (see Lemma [7. II below) . First, rewrite ( 14. ip with q = q^ as 

A D (x) = R p QD (x)/b D (n) - R VA D (x) , (5.2) 

where Ad(x) := Be>(x) /br>{n). Then the first term on the right-hand side of (15.21) converges to 
the first term on the right-hand side of (15.11) : 

Rop QD (x)/b D (n) -> e ikx , D ->■ oo , (5.3) 

in C(R 3 ) due to (@3D and (12131 . 

It remains to prove the convergence of the corresponding solutions: 

Proposition 5.1. Let conditions H1-H3 hold and k ^ 0. Then 

A D (x) ->■ , D — )■ oo , (5.4) 

&n "Hlo- wf/i any s < 2 and cr > 5/2. The function A(x) is defined by (II. lip . 

Proof. We deduce the convergence from the compactness of the family {Ad(x) : D > 0} and 
the Ikebe Uniqueness Theorem [31 Theorem 3.1], which holds under the condition k ^ 0. 

Step i) Formula (I5.2p implies 

H^lki < C\\Rop qD \\ H 2_ +\\R VB D \\ H 2_ . 
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The first term on the right-hand side is uniformly bounded for D > since estimates of type 
(14 . 1 Op — f 14 .lip hold with (x)~ u instead of V(x). The second term is uniformly bounded since VB q 
is uniformly bounded in C? a with a < 5/2 + e by (14. 7p while the operator R : C 2 a — > W?_ a is 
continuous for a G (1/2, 5/2 + e). Indeed, (— |A — E k )R Q = I, and hence the operator 

(_I A + 1)R — (I + E k + 1)R : Cl C 2 _ a 

is continuous by [El Lemma 2.1]. Hence, 

sup||A D || W 2 < oo , a G (5/2, 5/2 + e) . (5.5) 

D>0 

Step ii) Now the Sobolev Embedding Theorem [T7] implies that the family {Ae>(x) : D > 0} is 
a precompact set in the Hilbert space % s _ cr with any s < 2 and a G (5/2,5/2 + e). Hence, for 
any sequence Dj — > oo there is a subsequence L'j/ — >■ oo such that 

A Dj ,( x ) ->A<0) , j'->oo, (5.6) 

where the convergence holds in Ki^ with any s < 2 and a > 5/2. Therefore, 

^ .,(x)^U(i), j'^oo, (5.7) 

where the convergence holds in H s a with s < 2 and some a > 5/2 by H3. 

S'tep iii) At last, equation (15. 2j) and convergences (15. 6j) . (15. 7p . and (15.31) imply the equation ( 15. ip 
for 

A,(x) = e ifcx - J R ^(x) (5.8) 

since the operator i? := Ro(E k + i0) : £^ — >■ £ 2 CT is continuous for a > 1/2 by [161 Lemma 2.1]. 

The function A*(x) is bounded by (I4.2ip . and continuous by the Sobolev Embedding Theorem 
since A*(x) G W_ a with any s < 2 and a > 5/2 by (15. 6p . 

Finally, A(x) = A*(x) by the Ikebe Uniqueness Theorem [31 Theorem 3.1] for bounded 
continuous solutions to the Lippmann-Schwinger equation (15.11) under condition (12. 3p . Hence, 
convergence (I5.6P implies (15 .4p since the limit function A* (x) does not depend on the subsequence 
j>. □ 

Remark 5.2. Let us emphasize that the right-hand side of (I4.2ip with q = —nD is not uniformly 
bounded for D > : its value at x = q tends to infinity as D — > oo. 

6 Convergence of flux 

We check the convergence of the limit flux as the source goes to infinity. First we check (12. lip 
and (12TT21 relying on (EZJ) and ( 12101) . 

Lemma 6.1. Under conditions HI - H3 the long range convergence 

ip D (x,t)/b D (n) ->• A(x)e~ iEkt , D oo , ieR, (6.1) 

/io/ds in "Hi CT wit/i any s < 2 and a > 5/2. 
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Proof. The convergence follows from Proposition 15.11 since 

<Pd(x, t)/b D {n) = A D (x)e- lEkt (6.2) 

as t — > oo by definition of ip q (x,t) in (j2.7p with q(D) = —nD (recall that we omit the sum over 
the discrete spectrum in (]2.7]) ). Here \br){n)\ = \b(n)\ ^ by the Wiener condition (12. 4p . □ 

Lemma 6.2. i) Convergence (I2.12p holds: 

j D (x) = lm[[<p D (x,t)V(p D (x,t)] ^joo(x) = \b(n)\ 2 lm[A(x)VA(x)} (6.3) 
as D ->■ oo in Cf oc [R 3 ). 

ii) Moreover, the convergence holds "in the sense of flux", i.e., 

I 3d{x) • v{x)dS(x) I joo(x) ■ u(x)dS{x) , D oo , fel, (6.4) 
Js Js 

for any compact smooth two-dimensional submanifold S C M 3 with a boundary, where u(x) is 
the unit normal field to S, and dS(x) stands for the corresponding Lebesgue measure on S. 

Proof. Convergence (16.11) also holds in C(1R 3 ) since H s _ a C C(1R 3 ) with s > 3/2 by the Sobolev 
Embedding Theorem [17]. Further, the convergence of the derivatives 

Vtp D (x,t)/b D (n) ^VA(x)e~ iEkt , D ->■ oo , teR, (6.5) 

holds in H^T 1 with any s < 2. Hence, convergence (16. 3p holds in ZZj^IR 3 ), and moreover, 

V<f D {x,t)/b D (n) s ^VA(x)e~ iEkt s , D ->■ oo , teR, (6.6) 

in C 2 (S) by Sobolev Trace Theorem [T7J since we can take s — 1 > 1/2. 

Therefore, the integrands in (16.41) converge in C 2 (S) by (16. 3p . □ 

7 Long range asymptotics 

We obtain asymptotics (I1.12p . The first lemma is well-known |24j . 

Lemma 7.1. Equation (15. ip admits a unique bounded continuous solution which is given by 

dm]): 

A(x) = e ikx - RVe ikx (7.1) 
Proof. We should prove (17.11) assuming (15. ip . First, we apply the general operator identity 

P- 1 = Q- 1 + Q-\Q-P)P- 1 
to P = Hq — Ek — iO and Q = H — Ek — iO. Then we obtain Rq = R + RVRq, and hence 

RoVA = RVA + RVRqVA = RV(A + R VA) = RVe ikx 
by dSHJ. Substituting into ([51]), we obtain (173]) . □ 
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Next we need a generalization of Lemma H~T1 to functions of weighted Agmon-Sobolev spaces. 
Lemma 7.2. Let r(x) G TL 2 , where a < 7/2 + 5 with a 5 > 0. Then 

e i\k\\x\ 

R^rix) = (f)(6)— — + K(x) , xeR 3 , 9 := — , (7.2) 
\x\ \x\ 

where 

<t>{9) = Vtort(\k\9) , \9\ = l, (7.3) 

and 

\K(x)\ + \VK(x)\ < C\x\' 2 . (7.4) 
Proof. We adopt further the arguments from the proof of Lemma 3.2 in [3J p. 275]. First, 

e i\k\\x\ r x i r ly\2 

Ror(x) = ——- e M ] w' y r(y)dy + — / 1 -B(x,y)r(y)dy , 

2tt\x\ J \x\ J \x-y\ 

where the function B(x,y) is bounded. Hence, it remains to check that 



J(x) ■= j ^ \r(y)\dy < C\x\- X . 



\x - y\ 

Using the Cauchy-Schwarz inequality, we obtain 

( f 1 \ 1/2 

|J(x)| -ii R^P ^ l|r|U| - 

Now it suffices to prove the bound 

1 

b-j/Ry) 2 *- 4 

In the spherical coordinates, similarly to (14 . 1 6 j) . we obtain 

r 2 dr r sin 9d9 



Hz) ■= / , .,, 9 , A ^ dy < C\x\- 2 . (7.5) 



f°° r 2 dr f 
I(x) = 2tx I - --— / 



o (l + r) 2<T 4 J \x\ 2 + r 2 — 2\x\r cos 9 



o I I / pdp i \ l + P\ 

MA L (i + pN)— log ^ 



< c r , n ^_, + cixi 5 -- + c r wp 



o (l+p|x|) 2 - 4 1 1 J 3/2 (l + p\x\) 



n , > l lXl/2 r2dr ,5 2, „ dT 

c \ x \ / h , + ^M 5 ~ 2,J + c 



2ct-4 



(1 + r) 2 - 4 1 1 7 3 N/2(l + r) 2CT - 4 



< CUaMxl- 2 + C 2 \x l5 - 2a 
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if 2cr > 7. Hence, (I7.5P holds since we can take any a < 7/2 + 5. Then the first bound in (I7.4p 
follows. To prove the second bound, we differentiate (17. 2\i : 

i\k\\x\ 

VR o r(x) = 0(6)i\k\6^—— + O(\x\- 2 ) + VK(x), \x\ -> oo . (7.6) 
On the other hand, VRor(x) = RqVv(x), where Vr(x) G T-L\. Hence, by the arguments above, 

i\k\\x\ 

VR r(x) = 1 (0)i— - + 0(\x\- 2 ) , bl^oo, (7.7) 
\x\ 

where 

^(0) = V2^\Tr(\k\6) = V27ri\k\9r(\k\9) = i\k\9<j>(9) . 
Hence, the second bound in (I7.4p follows by comparison of (17. 6p and (17. 7p . □ 

Now asymptotics of type (11. 12ft follow from ( 17. ip and the next lemma, which is a modification 
of Theorem 3.2 from [3J Chapter 4]. 

Lemma 7.3. Let condition H3 hold. Then 

-RVe lkx (x) = a(k,9)^-^- + KAx) , Id ->■ oo , (7.8) 

\x\ 

where 9 := the amplitude a(k,9) is given by (I1.15p . and 

\Ki(x)\ < C\x\~ 2 , |V#i(z)| < C\x\- 2 . (7.9) 

Proof. First, RV = R T, where T := T(u + iO) (see (3.31) of Chapter 4] and 0); hence 

- [RVe lkx }{x) = -[R Te tkx }(x) . (7.10) 

We cannot apply Lemma 14.11 directly to obtain asymptotics (I7.8P for R Te lkx since we did not 
prove a bound of type flM} for Te ikx . Let us prove a weighted estimate for Te lkx and then apply 
Lemma 17.21 

Namely, Te ikx G H% with any a < 7/2 + e/2. Indeed, 

Te ikx = Ve ikx - VRVe ikx , 

where Ve ikx G H 2 a with a < 7/2 + e/2 by H3. Furthermore, V-RVe ifcB G ft* wit h a < g/ 2 + e 
since #e ifa G ft* wit h o- < -1/2. The latter follows from RVe lkx G C 2 a with a < -1/2 by the 
identity 

(1 - -A)R{u) = 1 + (1 - V + u)R{u) . 

Finally, applying Lemma I7T21 to the function r(x) = Te lkx and using (I7.10p . we obtain asymp- 
totics (JZSD with 

a(k, 9) = -V2^r(\k\9) = -(27r) 2 (Te fc , e w ) = -^ix 2 T{\k\9, k) (7.11) 
according to ( II. 7p . Hence, formula ( I1.15P is proved. □ 
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8 Cross section 

Now we can justify formula ( jl,14p . Convergence (16. II) can be written as the long range asymp- 
totics 

Lp D (x,t)~b D (n)A(x)e- iEkt = b D (n)e lkx e- lEkt -b D (n)RVe ikx e- lEkt , D > 1 , 

according to (17.11) . Hence, asymptotics (I7.8P and (I7.9P imply that the incident and scattered 
waves should be identified asymptotically as 

?p in (x,t) = b D (n)e ikx e- iEkt , ?p sc (x,t) = -b D {n)RV e ikx e~ iEkt , D > 1 . 

Then the corresponding limit stationary flux is given by 

j in (x) = \b(n)\ 2 k , f c (x) = \b(n)\ 2 lm[^{x)Va sc (x)] , 

where a sc (x) := — RVe lkx . At last, (I7.8P and (17. 9p imply that the scattered flux admits the long 
range asymptotics 

f c (x) = \b{n)\ 2 ^^^ \k\e + 0{\x\- 3 ) , \x\ -> oo , (8.1) 

with corresponding angular density (11.21) at infinity 

jl%9) = \b(n)\ 2 \a(k,x/\x\)\ 2 \k\. (8.2) 
Finally, according to definition (II. ip . the differential cross section is given by 

a(6):=j: c (e)/\r\ = \a(k,x/\x\)\ 2 , 
which justifies (11.141) . Then (11.61) also holds according to known formula (I1.15p . 
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